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Study of the hadronic production of kaon pairs below the threshold for 
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The near threshold production of K^K~ pairs in proton-proton collisions has been investigated at the cooler 
synchrotron COSY below and above the threshold for the (p meson using the COSY-11 and ANKE facilities, 
respectively. The excitation function determined for the pp — > ppK'^K~ reaction revealed a statistically significant 
enhancement close to the threshold which may plausibly be assigned to the influence of the K~p interaction. In 
addition, observed consistently by both groups, a strong enhancement at low values of the ratio of the K~p to K'^p 
invariant mass distributions shows that the proton interacts much stronger with K~ than with K^ . In this report 
we focus on the measurements performed by the COSY-II collaboration. We explain the experimental method 
used and present main results of completed analyses as well as a new qualitative elaboration of the ppK^K~ 
events population on the Goldhaber plot. We conclude with the observation that event densities increase at the 
region where the influence from the K^ K~ interaction is expected. 



1. Introduction 

A primary motivation for measuring cross sec- 
tions for the pp — > ppK^ K~ reaction near 
the kinematical threshold was the study of the 
hadronic interaction between K^ and K^ mesons 
in order to understand the structure of the scalar 
resonances /o(980) and ao(980) 1^. Such mea- 
surements have been made possible by beams of 
low emittance and small momentum spread avail- 
able at storage ring facilities and in particular at 
the cooler synchrotron COSY placed in Jiilich, 
Germany [2] . A precise determination of the col- 
lision energy, in the order of fractions of MeV, 
permitted to deal with the rapid growth of cross 
sections [3 and thus to take advantage of thresh- 
old kinematics like e.g. full space phase cover- 



age achievable with dipole magnetic spectrom- 
eters rather limited in geometrical acceptance. 
Early experiments on K^ K~ pair production at 
COSY conducted by the COSY-11 collaboration 
revealed, however, that the total cross section at 
threshold is by more than seven orders of magni- 
tude smaller than the total proton-proton produc- 
tion cross section making the study difficult due 
to low statistics |4|5|6j . A possible influence from 
the /o or ag mesons on the K^ K^ pair produc- 
tion appeared to be too weak to be distinguished 
from the direct production of these mesons based 
on the COSY-11 data [5]. Recent results obtained 
by the ANKE collaboration with much higher 
statistics can also be explained without the need 
of referring to the scalars /q or oq [718] . However, 
the systematic collection of data below |4I5I6| and 



above [719] the meson threshold combined to- 
gether reveal a significant signal in the shape of 
the excitation function in which the K~p and per- 
haps also the K^K~ interaction manifests itself. 
These observations motivate us to search for a sig- 
nal from the interaction between the K^ and K^ 
in two dimensional invariant mass distributions. 
The analysis is based on generalizations of the 
Dalitz plot for four particles proposed by Gold- 
haber et al. [10111] . The knowledge about the 
KK and KN interactions is important in many 
physical fields. In addition to the already men- 
tioned studies of the nature of the scalar resonaces 
ao(980) and /o(980), in particular for their inter- 
pretation as KK molecules |12I13I14] . it is also 
of importance in view of discussions on the struc- 
ture of the excited hypcron A(1405), since it is 
not clear whether it is a usual three quark system 
or whether it is a KN bound state p^ . Further- 
more, an understanding of kaon and antikaon in- 
teractions with a nucleon is essential for studies of 
properties of strange particles immersed in dense 
baryonic matter 16] and in the determination of 
the structure of neutron stars |17|18j . 

2. Measurements of the pp -^ ppK^K^ re- 
action at COSY- 11 

The measurements of the pp —> ppK^K~ reac- 
tion close to threshold have been conducted using 
the cooler synchrotron COSY [5] and the COSY- 
11 detection system il9i shown schematically in 
Fig. m The target, being a beam of H2 molecules 
grouped inside clusters of up to 10^ atoms |20) . 
crosses perpendicularly the beam of protons cir- 
culating in the ring. If at the intersection point 
of the cluster target and COSY beam a collision 
of protons leads to the production of a K'^K^ 
meson pair, then the reaction products having 
smaller momenta than the circulating beam are 
directed by the magnetic dipole field towards the 
COSY- 11 detection system and leave the vacuum 
chamber through a thin exit foils [19]. Tracks 
of positively charged particles, registered by drift 
chambers, are traced back through the magnetic 
field to the nominal interaction point leading to 
a momentum determination. Knowledge of the 
momentum combined with a simultaneous mea- 
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Figure 1. Schematic view of the COSY-11 de- 
tector with an exemplary event of the reaction 
channel pp -^ ppK^K~ . For the description see 
text. 



surement of the velocity, performed by means of 
scintillation detectors SI and S3, permits to iden- 
tify the registered particle and to determine its 
four momentum vector. Since at threshold the 
center-of-mass momenta of the produced parti- 
cles are small compared to the beam momentum, 
in the laboratory frame all ejectiles are moving 
with almost the same velocity. This means that 
the laboratory proton momenta are almost two 
times larger then the momenta of kaons. There- 
fore, in the dipole field protons experience a much 
larger Lorentz force than kaons. As a conse- 
quence, in case of the near threshold produc- 
tion, protons and kaons are registered in separate 
parts of the drift chambers. Therefore, as a first 
step in the reaction identification events with two 
protons registered in an appropriate part of the 
drift chambers are selected based on the time- 
of-fiight between the SI and S3 scintillation ho- 
doscopes. The additional requirement that the 
mass of the third particle, registered at the far 
side of the chamber with respect to the circulat- 
ing beam, corresponds to the mass of the kaon, 
allows to identify events with a. pp —t ppK'^X~ 
reaction [16] . Knowing both the four momenta of 
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Figure 2. (Upper panel) Exemplary missing mass 
spectrum determined for the pp -^ ppK^ X^ re- 
action at an excess energy of Q = 17 MeV [S]; 
(Lower panel) Missing mass distribution from up- 
per panel after additional requirement of the sig- 
nal in the dipole detector as it is described in the 
text. 



positively charged ejectiles and the proton beam 
momentum one can calculate the mass of an un- 
observed system X~ . Figure [5] (upper panel) 
presents an example of the missing mass spec- 
trum with respect to the identified ppK^ sub- 
system. In the case of the pp — > ppK^K~ re- 
action this should correspond to the mass of the 
K~ meson, and indeed a pronounced signal can 
be clearly seen at this position. The additional 
broad structure seen in the upper panel of Fig.[2]is 
partly due to the pp ~f ppir'^ X~ reaction, where 
the 7r+ was misidentified as a K^ meson, and 
in part due to the K^ meson production associ- 
ated with the hyperons A(1405) or S(1385) [5IS] . 
The background, however, can be completely re- 
duced by demanding a signal in the silicon pad 
detectors (mounted inside the dipole) at the po- 
sition where the K~ meson originating from the 
pp — > ppK^ K^ reaction is expected (lower panel 
of Figure 2). This clear identification allows to se- 
lect events originating from the pp -^ ppK~^K~ 
reaction and to determine the total and differen- 



Figure 3. Total cross section as a function of the 
exess energy Q for the pp — *■ ppK^K" reaction. 
The data are form references [4I5I6I7I9] and the 
meaning of the lines is described in the text. 



tial cross sections. 

Figure [3] shows the excitation function for the 
pp — > ppK^K~ reaction established by the 
COSY-11 [4I5I6J group near the threshold and 
by ANKE [7] and DISTO [9] collaborations at 
higher energies. The dashed line shows the re- 
sult of calculations under the assumption of a ho- 
mogeneous phase space population, normalized 
to the DISTO data point (Q = 114 MeV). For 
the dashed-dotted line, the proton-proton final 
state interaction is included using parameteriza- 
tion known from the three body final state J21] . It 
is clearly seen that calculations neglecting the in- 
teraction of kaons underestimate the experimen- 
tal results by a factor of five in the vicinity of the 
kinematical threshold. Therefore, the enhance- 
ment may be due to the influence of K^p or 
K^K~ interaction. And indeed, with a factor- 
ization ansatz for the pp and pK^ interaction, 
the ANKE collaboration described the excitation 
function much better [718] , however still underes- 
timating the two lowest data points by more than 
a factor of two. This could indicate that in this 
region the influence of the K^K~ interaction is 
significant and cannot be neglected. This obser- 
vation encouraged us to carry out the analysis of 



differencial cross sections for tlie low energy data 
at Q = 10 MeV (27 events) and Q = 28 MeV 
(30 events) , in spite of the fact that the available 
statistics is quite low [B]. 

3. Goldhaber plot analysis: generalization 
of the Dalitz plot for four particle final 
state 

Usage of the Dalitz plot for extracting informa- 
tion about the interaction among particles in the 
case of three body final states is well known. It 
was introduced by Dalitz in a nonrelativistic ap- 
plication |22j and then extended to the relativis- 
tic case by Fabri [23] . If the transition amplitude 
is constant over phase space and if additionally 
there is no final state interaction, the occupation 
of the Dalitz plot would be fully homogeneous 
because the creation in any phase space interval 
would be equally probable. Thus, final state in- 
teraction should show up as a modification of the 
event density in the Dalitz plot. 

In the case of four particles in final state the 
analysis is more complex, because one needs five 
variables to fully describe a relative movement of 
particles. Nevertheless, there are many different 
types of generalization of the Dalitz plot for four- 
body final states. In this contribution we present 
a generalization proposed by Goldhaber [10|llj . 
which we use further on for studying the inter- 
action in the ppK^K^ system. However, there 
exist many other approaches as described e.g. by 
Nyborg or Chodrow [24125126] . Consider a reac- 
tion yielding in the final state four particles with 
masses rrii and total energy y^ in the centre- 
of-mass frame. Assuming that the matrix ele- 
ment for the process M depends only on invariant 
masses of two- and three particle subsystems [23] 
the distribution of events can be expressed in 
some choice of five independent invariant masses 
e.g.: Mf^, M|4, Mf^, Mf^^, M^^^. Assuming 
further that M depends at most on M12, -^34, 
and M^24j which corresponds to the situation 
where only two two-particle or one three-particle 
resonances are present [24], one obtains the fol- 
lowing distribution of events: 
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Figure 4. Goldhaber plots for the pp -^ 
ppK^ K" reaction, (upper panel) Simulations at 
Q = 28 MeV for homogeneously populated phase 
space modified by proton-proton final state inter- 
action taken into account as weights proportional 
to the inverse of a squared Jost-function of the 
Bonn potential [21127] . (middle and lower pan- 
els) Experimental data obtained at Q = 10 MeV 
(middle), and Q = 28 MeV (lower). 
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where g is a simple analytical function of its vari- 
ables [21j. The projection of the physical region 
on the {Mi2, M34)-plane gives a right isosceles 
triangle in which the area is not proportional to 
the phase space volume. It is important to note 
that the event density in the Goldhaber plot is 
not homogeneous and goes to zero on the entire 
boundary of the plot given by the following equa- 
tions: M12 + M34 — y^, M12 — mi +1712, M34 — 
1713 + 1714 ED- 

Figure |4] (upper panel) presents the simulated 
distribution for the ppK^K^ reaction deter- 
mined taking into account only the pp~FSI. Ex- 
perimental event distributions after acceptance 
corrections for both studied excess energies are 
shown in the middle and lower panels. Clearly 
the event densities in the experimental Goldhaber 
plots differ from the simulated spectrum. In par- 
ticular data show an enhancement in the range 
of small K'^ K~ invariant masses which may sig- 
nify a signal from the kaon-antikaon interaction. 
A similar enhancement is seen by the ANKE 
group [7] below the <j) meson mass in the K~^ K^ 
invariant mass distributions. As a next step in 
the analysis we will compare experimental data to 
the results of Monte Carlo simulations generated 
with various parameters of the K~^K~ interaction 
taking into account the pK-FSl, as described in 
references I7I8I. 
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